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Nonlocal quantum corrections to gravity have been recently proposed as a possible solution to the 
cosmological fine-tuning problems. We study the dynamics of a class of nonlocal actions denned by 
a function of the inverse d'Alembertian of the Ricci scalar. Power-law and exponential functions are 
considered in detail, but we also show a method to reconstruct a nonlocal correction that generates a 
given background expansion. We find that even the simplest terms can, while involving only Planck 
scale constants, drive the late time acceleration without changing early cosmology. This leads to a 
sudden future singularity, which however may be avoided by regularizing the d'Alembertian. We 
also consider the Einstein frame versions of these models. 
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£ ! i- INTRODUCTION 

The observed acceleration of the universe requires the presence of a tiny energy density that corresponds to a 
\Q ' cosmological constant of the order A ~ 10 -120 in Planck units. This has motivated the consideration of dynamical 
alternatives to the cosmological constant These alternatives usually introduce new fields Q, by either adding 
quintessence in the matter sector or by modifying gravity. In the quintessence models the constant term is promoted 
to a potential function A — ► m v (p 2 of the cosmon field tp [3|, and in the modified gravity models one instead considers 



a 



a curvature-dependent effective cosmological term R + A — > R+ /i 2 (™+ 1 ) / ' R n [J|. Typically one still requires an 
explanation for the tiny constants m^/i ~ 1CP 33 eV within both of these approaches (in addition to the vanishing 
of the usual A). Recently a new, nonlocal variant was suggested [5| in the form R + A — > R + Rf(0 where the 
| function of the inverse d'Alembertian acting on the curvature could involve constants of (roughly) the order 

^ . one in Planck units. This seems a worthwhile starting point. 

' Quantum loop effects are well known to generate nonlocal corrections to the effective action of gravity @ . Nonlo- 
&\ , cality arises naturally in string theory, and is considered as the most likely resolution of the black hole information 
paradox Q ■ Cosmological effects have been shown to arise from nonlocal corrections which were proposed as a cure to 
the unboundedness of the Euclidean gravity action Q. The resulting curvature dependence of the effective Newton's 
£C) • constant may then easily pass the local gravity experiments, but could be constrained by its cosmological effects, 
especially on nucleosynthesis. String-inspired forms of nonperturbative gravity could also be asymptotically free and 
thereby describe bouncing cosmologies [1, [l(| • The cosmological impact of the particular (logarithmic) form of quan- 
tum corrections expected at the leading order has been computed and a small possible effect on the expansion rate 
, was found (lH Il2| . However, since it is a nonrenormalizable theory, long-distance modifications of gravity which could 
become significant at low curvature cannot be excluded, and in particular they could be responsible for the observed 
rN | acceleration of the universe [3, [3] ■ 

Let us note also that nonlocal cosmology has gained recently interest in the context of string field theory, with 
various applications to incorporate bounces, model dark energy and generate nongaussianities during inflation Il5l . 
[H, [13, Us, [H, H3, HH IHl ■ For recent clarifications about the degrees of freedom in nonlocal field theories, see [H, |24| . 

In this paper we focus on the cosmological consequences of the class of actions describing additions to general 
relativity as a function / of the inverse d'Alembertian (i.e. covariant Laplacian) acting on the curvature (i.e. Ricci) 
scalar. The dynamics of the cosmological background expansion is studied with the main aim to find solutions 
naturally incorporating the present acceleration of the universe. Then it is crucial to uncover in detail whether these 
models could feature a viable sequence of radiation dominated, matter dominated and accelerating eras. In particular, 
we then study the dynamics in two simple cases: when the function / is of the power-law and of the exponential form. 
We are then able to find the conditions for the late acceleration to occur, and to determine the parameters of the model 
required for observationally consistent evolution of the universe. We also consider an alternative approach, where we 
decide on the background cosmological expansion history and reconstruct the form of the function / to generate the 
given expansion. In general, such a procedure is difficult to perform analytically, and it may yield functions which are 
very complicated, unstable or singular. Therefore we mainly consider the approach of choosing a simple form of the 
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function and studying its consequences, which can be then uniquely determined. In fact, we find that the evolution 
is qualitatively similar regardless of whether / is exponential or a power-law. 

In these models the acceleration of universe does not lead into a de Sitter stage, but instead ends in a singularity, 
where the pressure diverges at a finite time, when the scale factor has expanded a finite amount. This is a so called 
future sudden singularity [25| . which is different from the big rip singularities present in phantom dark energy models, 
described by a diverging scale factor at some finite moment of time. According to the classification of Ref.[26|], the 
sudden future singularity is labeled type II. We show that in the singularities may be avoided by regularizing the 
d'Alembertian operator. 

The plan of the paper is as follows. In Section [TT] we write down the basic equations, introducing the model and 
it's presentation as a biscalar-tensor theory in III Al and writing down the cosmological equations in III Bl We also 
discuss the Einstein frame version of the theory in III CI In Section IIIII we then study cosmological dynamics from 
various points of view. The initial conditions are conceptually important in these theories, and in this particular 
class of models we can set the initial conditions at the radiation-dominated era according to the analytical solutions 
presented in subsection IIII Al We find also analytically the conditions for the matter-dominated era to end. These are 
derived in subsection IIII Bl for the power-law and the exponential model. Numerical computations show the details 
of the evolution and reveal the generality of the singularity. The approach to singularity and its avoidance is studied 
analytically in a particular case in subsection IIII CI An alternative route to finding solutions is given in subsection 
IIII Dl as a reconstruction method and it is applied in a discussion of scaling solutions in the special case of exponential 
/. In Section ITVl we present our conclusions. 



II. THE NONLOCAL MODEL 

In this Section we specify our nonlocal action and rewrite it as a multi-scalar-tensor theory [2tJ with extra kinetic 
couplings. The cosmological equations are then written as generalizations of the standard Friedmann and continuity 
equations, and this is again recasted into more useful form by a change of variables. Finally, the corresponding system 
is considered in the Einstein frame. 



A. Multi-scalar-tensor theory 

We consider the following simple example of a nonlocal action [28| 

S = J d n xj—g [R (1 + /(D^i?)) + 2n 2 C m ] . (1) 

Introducing the Lagrange multiplier £ to rename the inverse dAlembertian as 0, we may rewrite the action in the 
form 

S = J (FxV^g [R (1 + /(</>)) + £(□</> -R) + 2n 2 C m ] . (2) 

More general action involving m different powers □ acting on R can be recasted into a multi-scalar-tensor theory by 
introducing 2m scalar fields [28j |. Here we restrict to the case where only the first power D _1 appears in the action. 
Integrating partially and neglecting the boundary terms, as usual, we can then eliminate the double derivative, 

S = J d n x^—g [R (1 + /(</>) - - V Q £V Q + 2k, 2 C m ] ■ (3) 

This is a peculiar kind of double-scalar-tensor theory, where massless, nonminimally coupled scalars have a kinetic- 
type interaction resembling of the nonlinear sigma model. From the form ([3]) it is clear that one may not eliminate 
the kinetic coupling by a field redefinition. However, it turns out to be convenient to make yet one more manipulation 
of the action. If we use, instead of the field £, the field 

o = m -e (4) 

as our dynamical variable, it becomes transparent that we may consider a double-scalar-tensor theory, with only one 
of the scalar degrees of freedom is nonminimally coupled to gravity. The pair of the scalar fields still has the unusual 
kinetic- type coupling [i^. 

S = J d n x^—g [(1 + *)R - f{<j>){d(j>f + V Q *V> + 2n 2 C m ] . (5) 
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We get the field equations 

Gftv = s^-iVaCV'^ + V^V^-^n-V^-G^JC/C^-O (6) 
= tJT^+T^+T^+T*, (7) 
where the effective energy momentum tensors for the two fields and for their coupling read as 

1 



T+ v = /'(</>) 



(8) 



T* = -(G^ + g^n-V^Vv)*, (9) 

T%, = ^aW^-V^)*. (10) 

The first piece reduces to canonical kinetic terms, when f = <fi. Because of the additional coupling to the other 
field, /'(</>) > does not now necessarily imply that the field is a ghost. The second part corresponds to the usual 
nonminimal coupling of a standard Brans-Dicke field \&, and is also well-known to be consistent contribution to field 
equations both at the quantum and classical level given W > — 1. The last piece is linear in the derivatives of both 
the scalar fields, which may not be an unproblematical feature. 
The equations of motion are 

□* - /"(0)W) 2 +2/'(0)i?, (11) 
□ = R, (12) 

where the Ricci curvature could be solved from the trace of field equations, 

(1 - |)(1 + *)R = k 2 T - (2 - l) (/'(0)(^) 2 - V Q W^) - (»- 1)D*. (13) 

Since minimally coupled, the matter tensor is conserved, V M T™, which may be used to check the system of equations 
(f7l TT3|) as it is not an independent condition (29|. 

Let us consider the case / = foR/U. From the equations of motions (fTTj) one sees that U<p is now proportional 
to the It is thus possible to eliminate one degree of freedom. This is similar to the so called nonlinear gravity 
theories defined by a function of R: in such a class of theories there is an extra scalar degree of freedom, which vanishes 
only in the unique case of the Einstein- Hilbert action that is linear in R. Here we then noticed that multiplying the 
curvature term by a function of OR is equivalent to introducing two additional scalar fields, and the unique case of a 
linear function again allows one to integrate out one scalar degree of freedom. 

One could also consider generalizing the D'Alembertian operator □ with regulating terms like □ — eR + 
(~f/n 2 )R^ prT R flup ' J ■ ■ ., where e and 7 are dimensionless numbers [8(. If we make the substitution □ 1 — > (□ — h(R)) 1 
in the original action |T]), by similar steps as before we end up with a generalized Lagrangian, 

d"aV=S [R(l + /(</,) - £ - h'( X m + (h'(x)X ~ h( X m<t> - V Q £V> + 2 « 2 £ m ] . (14) 

If one considers a linear function h{x), it is possible to eliminate the degree of freedom x- Even with a general h(x), 
in the case of linear /(</>), one may eliminate one degree of freedom since £ can be related to <fi. Thus, with linear 
/ and linear h, Eq. (fT4]) becomes an action for a single scalar nonminimally coupled to gravity. We return to this 
interesting special case in MI CI 



B. Cosmological equations 

Having now the general equations at hand, we may see how they become in cosmology. In a flat Friedmann- 
Robertson- Walker (FRW) background, 

ds 2 = -dt 2 + a 2 (t)dx 2 (15) 

we get the Friedmann equations 

3H 2 (1 + tf) = K 2 Prn + X - (/'(^)0 2 - iffy - 3Hi> (16) 
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- (2H + 3H 2 ){1 + *) = K 2 Pm + ~ (/'(0)0 2 - tf^) + * + 2HV (17) 

The equations of motion for the fields are then 

# + 3i3> = -6(H + 2H 2 ), (18) 

* + 3iT* = f"(<j))<j> 2 - V2f{(j)){H + 2H 2 ). (19) 

Introduce the dimensionless variable 

together with the e-fold derivative y* = dy/d\oga. We define the energy density fraction of matter as usual, 



O m ee ^L, (21) 



We then have a complete set of dimensionless variables. The dynamical system may then be written as 

fl* n = 3(w m -w e ff)Q m , (22) 
*** = ^K//-l)^-6/'(0)(l-3 We// )+4r(^)X 2 , (23) 

X* = ^[(w eff -l)X + 3w eff -l}. (24) 

where 

3w m n m + 2/' (0)(X 2 - 3) + 4/"(0)X 2 - (1 + X)^* 
= 3(l + *-6/'(fl) • (25) 

The first equation simply expresses the usual matter conservation law, and the two following equations are the 
conservation laws for our two scalar fields. 

The fixed points of the system, defined by fl* n = = X* = are as follows. 

• It is possible to construct scaling solutions with w e f / = w m and with arbitrary f2 m if the coupling is of the form 
/(</>) ~ e Ql *, where a = 3(1 — w m ) 2 / (2 — Qw m ). However, this is a very special case with not only the form but 
the parameters of the coupling fixed. Therefore we do not consider this possibility further. 

• In the case of effectively complete radiation domination, w e ff = 1/3, we automatically have a fixed point, since 
with X = the equations are identically satisfied. The matter density and the field $ are then related by the 

= 1 + In the very early universe, when we have exactly w e ft = w m = 1/3, the correct initial condition 
is VP = 0. The evolution begins when w m drops below 1/3, as we discuss in the following subsection. 

• If /' = /" = the equations are again identically satisfied with arbitrary 'J, and now X = 3(1 — w e //) 2 /(2 — 
6w e ff) and 1 + = Q m , The most relevant case is of course ^ = 0, which corresponds to standard cosmology. 
The field cj> keeps growing, as we see in the following, but it depends on the form of the coupling whether this 
means that the effect of nonlocal gravity is becoming significant. 

In fact we note that the information about the fixed points in this system throws light only on some special corners 
of the phase space, and for more satisfactory understanding of the nonlocal dynamics we now turn into other consid- 
erations. The system of alternative variables normalized by f(<fi), say £l m / f{4>), would feature additional fixed points 
in the special case of exponential potential (3(| • We will also find and generalize these same solutions in subsection 

dnTDli . 

C. The Einstein frame 

Though we will work mainly in the usual Jordan frame, we will briefly consider the Einstein frame version of these 
models in this subsection. Perform the Weyl rescaling [3l], HH 

~ 9iiV = (l + *)^^ = e^% v . (26) 
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The action becomes 

S = I d n x^-g 



R- 



n — 1 , ~ 



n- 2 



(V<^) 2 - e- v /'(0)(V0) 2 + V Q <^V Q </) + 2 K 2 e-^£ m (5e-^). 



(27) 



Now has also acquired a kinetic term. The curvature coupling has been removed in this frame, but the matter sector 
has now become nonminimally coupled. Since there remains the derivative coupling term, things have not simplified 
so much as in the Einstein frame of the Brans-Dicke class of theories. However, the metric rescaling provides also 
here a useful alternative viewpoint to the system. The field equation is 



where we have again decomposed the left hand side to four contributions, 

1 ,^ 



f+ = f'We-* 



T = \ 1 



1 



n- 2 



(28) 

(29) 
(30) 
(31) 



Note that none of the four energy momentum tensors is conserved. Matter, as well as the field <p, is coupled nonmin- 
imally to the field (p. The equation of motion for this field may be written as 



/ e v 2(n-l) 
\JrW)~ (n-2) 



□ 



(V0) 2 -V Q ^V^ 



k T 



(f - l) e 2¥>/(«-2) ' 



(32) 



This could be derived either by varying the action (|27[) or by conformally transforming the equation (jTTJ) and using 



R = 



n-l,~ 



1 _ 2 

x 2 



T + ^(V^) 2 - e-*7'(^)(V0) 2 - e-vVfrW^ip. 



(33) 



Then, to find the equation of motion for the field <fi in the Einstein frame, it is most convenient to note from Eqs. (|lip 
that □</>=§(□*- /"(0)(V0) 2 )//'(0). 

The Einstein frame cosmology can then be obtained from the Jordan frame solutions by carefully taking into account 
the conformal dimension of each quantity. From the scaling (|26[) we see that the scale factor, being a squareroot of a 
component of the metric, has conformal dimension half, The invariance of the line-element then requires 

dt = e v ^ 2 dt. From this already follows that the Hubble rate in the Einstein frame is 



H = e-*l\H+~<p), 



and that its derivative may be written as 



1 dH 



1 



H 2 dt l + 2<p* 



2 + ip* 



H* 1 
IT 



(34) 



(35) 



This dimensionless quantity gives directly the effective equation of state in the Einstein frame. As a manifestation of 
the matter nonconservation, we have p m = e~ 2cp p. m . It follows that 



fi, 



ev(l + i</?*) 2 



(36) 



gives the relative amount of matter. 



III. COSMOLOGICAL DYNAMICS 



In this Section we will study cosmology at four different epochs: radiation dominated, matter dominated, acceler- 
ating, and a general scaling. We can then find different approximations which are helpful in analytical considerations 
in these four cases. Numerical solutions are given in table [J and the figures [1] and [3] 
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A. Radiation domination: Initial conditions 



Consider the radiation dominated epoch. We assume that the nonlocal corrections are negligible. By definitions, 
we may always write 

R = 3(l-3w eff )H 2 . (37) 
Assuming standard gravity with matter and radiation as sources, we have 

1 / H \ 2 tt° 

to, eff = -—, hr =-f(l + ro), (38) 
1 + ra \Hq J a 4 

where the fraction of nowadays relative matter density, fi^ = 3p°„/ (8ttGHq) and nowadays relative radiation density 
0° = 3p®/(8ttGHq) is denoted r = Q^/Q,^,. When all matter is relativistic, R vanishes identically. Thus 4> should 
also vanish. When part of the matter becomes nonrelativistic, R begins to evolve (even if the universe is dominated 
by radiation). The scalar field (j> is given as the solution to the equation <f> = OR. Explicitly, we get 

m ~ ~ II J® £ « a ( t ") R ('")<«" + *'/'(«Fl) S dt ' + *>■ < 39 > 

Note that this does not depend on the form of f{4>), since the scalar (j> is just a way of reexpressing R/D. We consider 
exclusively the retarded solutions to respect causality. In our models, we then set (f>i = cf>i = at an initial, completely 
radiation dominated time ij. When R vanishes, <p = OR should also together with its derivatives. Using the expansion 
(f38|) . we can then reduce the double integral (f39|) into the form 

da' 

a ,2 VT+ra' 7 Qs vT 



^(°) = - 3r / ,2 „ . , I n , „ ( 4 0) 



which one may perform analytically, yielding 



4(a - Oj) a — / (l + Vl + m)(l- y/l + rgj) 2^1 + ra 2y/l + r§ \ 

0(a) = 2 log — + (rai - 2 VI + ra. ( log + . (41) 

raai a t \ (1 - y/l + ra)(l + V 1 + rai) ra r o-i J 

Obviously, the solution depends on the moment aj at which we assume the matter in the universe to become nonrel- 
ativistic. One then obtains the derivative of the field, 



X(a)=[2 + {rat-2)^^^]—-l. (42) 



y/1 + rai J ra 

Using the solution for <fi, the other scalar field £ may then be computed as 

m J* a3 (nR(t'vm"))dt" +£< £ (43) 

The value of £ thus depends on the form of /(</>), and has to be therefore considered separately for each form of 
the nonlocal term /. Analytic solutions may be difficult to find, unlike for <f>. Once £ is then solved, it is however 
straightforward to obtain $ = / — and its time first derivative. These then completely fix the initial conditions for 
our system. 

In principle, due to the nonlocal property, the initial conditions could carry information from inflation or even from 
a pre-Planck era, if there effectively R ^ 0. Here we just assume that (f> begins to evolve from zero when matter 
becomes nonrelativistic. Furthermore, we have assumed above that this happens instantaneously at a given epoch 
a = cti (in our numerical examples this is a, = e -23 ). In reality each particle species begins to contribute to the 
trace T m when the temperature of the universe has dropped below the mass of the particle. Thus a more precise 
calculation should consider r replaced by r*(a) in Eq. (|38p and elsewhere, r*(a) being a function which is zero as a — > 0, 
but is "kicked" to a slightly larger value each time some particle becomes nonrelativistic, and thus only eventually 
reaches r*(a) — ► r. Here we have approximated this a function as the step function r*(a) = r6(a — Oi), since the 
thermodynamic details of the very early universe seem not to be important for our conclusions. This is because the 
solution (|42p is an attractor. One may note this easily by inserting X + SX to the Eg. ([24]) . which then yields 



TO* = ~(3--i-KY. (44) 
2 V 1 + ar 
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Thus the evolution of <j> tends to a similar track with the solution (|41[). A different initial condition will only set a 
different normalization for the field. However, since values of the field in Eq. (|4"Tj) will be of the order of <j> ~ 10 -7 in 
the early universe, these differences are negligible when the correction terms become important. Note however, that 
in the case where the nonlocal corrections would be large at high i?, one should be careful with the early radiation 
history. 



B. Matter domination: Conditions for growth of nonlocalities 

At late enough matter era, when a>fli, the field (|41 [) is approximated by 

#o) = -2(log^ + l). (45) 

This limit also confirms our claim that at late enough era the dependence on the initial conditions is negligible, since 
the di is canceled out from the limit (|45|) . The field derivative during matter domination follows immediately, 

X(a) = -I. (46) 

This means that during an c-fold n _1 _R drops by —2. To analyze the evolution of the other field, we must then specify 
the function /. We consider two cases: an exponential and a power-law form. 



1. Exponential model 



Consider the exponential case, f(4>) = fo& ■ One may note that the pure exponentials models are ruled out 
because they feature a constant shift of the effective Newton's constant by /o, since the constant (f + fo) / (2k 2 ) then 
multiplies the standard Einstein-Hilbert term. This however may be avoided by simple reformulations of the model, 
and the exponential case is still an interesting example. We can now write 



f(a) = foa 



fo=[- 

re 



2n 



(47) 



One can then easily see that the leading contribution in the integral 
only when a < The solution is then 



6 - 4a 
3 -4a 



-2a 



43]) is growing with a proportional to £ 



-2<v 



(48) 



This solution satisfies the evolution Eq.(|23|) in the limit w e ff = 0. Since the coupling is growing, the correction 
term will inevitably become significant at some point. However, whether the nonlocal energy density then becomes 
dominant with respect to dust density, depends on the model. The effective equation of state of the universe is 
approximated by 



u ef f « 4a(a - 1)/ O o 



(49) 



Thus, only if fo < 0, do we expect the nonlocal energy density to dominate the matter density. When this happens, 
our solution (|45p is no longer valid. The approximations for the two fields are compared to numerical results in the 
right panel of figure [TJ When a < 0, the numerical example for an fo > case is depicted in the right panel of figure 
[21 and numerical examples for /o < cases are depicted in the right panel of figure [3] 



2. Power-law model. 



Consider the power-law case, f((f>) = fo4> n - F° r convenience we approximate log jra+1 P3 log ra and vT-f 
Then the solution to the double integral (|43| with the lower boundary terms dropped may be written as 



C(a) = -fo 



2 log- 

ra 



3 (ra) '■ 



(50) 
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From the fact that the incomplete Gamma function behaves asymptotically as T(n, x) m e x x n x (l + O(-)), we see 
that the second term in Eg. (|5D|) scales as ~ (log ra) n ~ 1 . Therefore the first term, which scales as ~ (log ra) n , will 
eventually dominate. Taking into account the contribution from f(4>) which is equal to £'s leading contribution, the 
coupling \& reads 

,, ^ !/ '( 2l ^)"( 1 -4)' < 51 ' 

One notes that this satisfies the evolution equation Eq.(|23p consistently up to the order ~ (log—)" 2 . Again we 
have a very simple condition for the nonlocal contribution to grow during an Einstein-de Sitter phase: if n > 0, the 
effect of the nonlocal terms increases during matter domination, otherwise not. Yet one may check what happens 
when these terms cannot be compleletely neglected. The effective equation of state of our solution is approximated 
by 

Weff * (-l)"/o— n ( log M) - (l + . (52) 

If (— l)"/o > 0, we expect that the standard matter domination may be interrupted by the growing nonlocal contri- 
bution, but since dust dilutes more slowly than the effective nonlocal density, acceleration will not follow. The case 
(— l) n /o < finally represents a scenario where f2 m — > asymptotically, and the nonlocal contribution will eventually 
dominate the energy density. The approximations for the scalar fields are compared to numerical results in the left 
panel of figure [TJ When n > 0, the numerical example for a case with (— l) n+1 fo > is depicted in the left panel of 
figure [21 and numerical examples for the opposite cases are depicted in the left panel of figure [3J 



C. Dark energy domination: A sudden singularity and its avoidance 



Given the conditions above, the universe will enter into an accelerating era. From numerical experimentation, we 
find that in the above simple cases this always leads to a so called sudden future singularity [25|]. This is type II 
singularity in the classification of Ref.[26[. There the second derivative of the scale factor diverges at a finite time 
and when the universe has expanded a finite amount since a^. In addition the first derivative remains finite. In terms 
of effective fluid quantities this means that pressure goes to infinity, while the energy density is finite. 

To avoid evolution into such a singularity, several effects have been proposed in frameworks of other models. 
The backreaction of conformal quantum fields might appear near the singularity and ch ang e the evolution [33| . 
Nonperturbative quantum geometric effects can also resolve the future sudden singularity [34j |. Here we show also 
that by slight modification of the nonlocal model by allowing a regulator term h(R) the singularity may cease to exist. 

For definitiveness and simplicity, let us consider the special case of linear coupling and a linear regulator term, 

f(x) = —tx, h(x) — ex. (53) 

As discussed in the end of section ([II Bl) . the action (Q3J) allows to integrate out two scalars. Then ^ = — r<p. The 
cosmological equations now simplify to § 



3# 2 (I - rq 



1 

2^ 



= -~j? + 3(T-e4)Hi>, 



-(2H + 3H 2 ){l-T(^+^ecj) 2 ) = (i + e)</> 2 + (r - e 

4> + 3H<p = 6(t - e<p){2H 2 + H) 2 . 



(54) 

(55) 
(56) 



The matter sources have been dropped, which is at least near the singularity is a good approximation since there 
matter becomes negligible, f2 m — ► 0. We can write the Friedmann equation (|55[) as a first order differential equation, 



da 
a 



2 1-r. 

This can be integrated. In the case e = 0, the solution is 




l ( l_ r ^ + l £ 02) 

3 (r - e0) 2 



o_ = / (Tfo-l)(l + y)(l-yo) 
a Q V N>-l)(l + y )(I-y) ' 



,vo-y 



(57) 



(58) 



9 



Because of difficulties in inverting this relation, we proceed by considering cosmology as a function of the scalar field 
4>. From Eq. ([57|) , we can relate the Hubble rate to the time derivative of the field for each value of the scalar. Then, 
we can relate effective density p$ = 3H 2 to the square of the field derivative, given the value of the field. To find 
a similar expression for the effective pressure p e ff — — (21? + 3-ff 2 ), we must first use equation (|55|) and solve 4> by 
using (|56p. The full result of this algebra, taking into account a nonzero e, is 



-6(1 + + 5r)) + c(r, </>, e)e - 4^6(1 - 6c) (t - e<t>) 2 J6 



2</> , -2+er 

W<p = 3(-2 + 2r(0 + 3r) - + 12t) + 6e 2 ^ 2 ) ' (59) 
where 

c(r, 0, e) = -3(-8 + (1 + 2e(-7 + 24e))</> 2 + 4(7 - 24er</> + 48r 2 )). (60) 

Consider first the case that e = 0. We then see that the equation of state diverges at finite (f> cr it — 3r — p The 
relation (|5"8|) shows that this happens at a finite value of the scale factor, and the Friedmann equation (|55[) shows that 
the density is also finite at 4> = 4> C rit- Thus we have a future sudden singularity. 

We then restore the regulator term e > into the operator to show that then the singularity may disappear. Now, 
the denominator of (15911 vanishes if 



r~6er± V(6e~l)(2e-r 2 ) 

(Pent = T"2 ( 61 ) 

e — be 

One may adjust e in such a way that the denominator in (|59[) is always nonzero, i.e. that the equation for the critical 
4> cr it as given by equation (|6ip does not have real solutions. To achieve this when the absolute value of r exceeds 
|t| > l/\/3, one has to set 1/6 < e < t 2 /2. When the absolute value of r is smaller, |t| > l/v3, one has to set 
t 2 /2 < e < 1/6. Thus we have shown that there are values of e which completely remove the singularity. 



D. Reconstruction of cosmology 



One might approach the dynamical problem from the other way and ask what kind of modifications / would be 
needed to generate a given expansion history. Previously such reconstructions have been considered for scalar-tensor 
theories [351 . vector- tensor theories [3(1 |37| . coupled scalar theories [38[, f(R) theories [1§| and the Gauss-bonnet 
theory [4Q|; for reviews, see fillip. 

Since our action contains an arbitrary function /, one might think that by adjusting the function suitably one could 
model an arbitrary background expansion. The recipe for finding the function / can be given quite simply. A given 
expansion history a(t) is equivalent to a given function w e ff{a)- Therefore, the variable X(a) and thus the field 4>(a) 
can be found by integrating the equation of motion (|24p 

X* = I [(w eff - 1)X + 3w eff - 1] . (62) 

It is possible to derive a similar equation for the other field * as well, 

*** + ~(7 - 3w eff )V* + 3(1 - w eff )9 = -3(1 - w eff ) + 3a»(l - w m ). (63) 

Therefore, a given background expansion indeed determines uniquely the form of both the fields <p and ^ . One can 
write the Friedmann equation in the form 

l + ^ = Q m + ^{2f'{ ( f>)X-^*)X-^*. (64) 

It is then straightforward to insert the solutions for cf> from Eq. (f6"2")) and the solutions for 'J from Eq. (f6"3")) into the 
Friedmann Eg. (f6"4")l to find f'{4>) as a function of the scale factor a. The coupling as a function of the scale factor is 
then given by the integral 



f(a)= / 2aX(a)f'(4>)da. 

J a,i 



(65) 
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/ Ww) = — o,.A 2 2 ( 68 ) 



If one is then able to invert 4>(a) to get a((f>), the functional form of the coupling is finally found as f(4>) = f(a(4>)). 
We demonstrate this procedure with an explicit example below. 

Consider cosmological expansion with the constant effective equation of state u> e // = w. We can then easily solve 
4> from Eq.dH]), 

0(o) = 2(1^) + + (g6) 

w — 1 

and check that the result also agrees with the integral (|39p . 0i and 4>o being the integration constants. In equation 
(|63| we set now f2 m (l — io m ) = fl^ n a 3 ^ w ~ Wm ^ (1 — ui m ). The equation is then solved by 

*(«> = -i + (^K««<"~ > + S + W 

where the two first terms give the special solution, simply a constant, and the two last terms give the general solution 
to the homogeneous equation with the constants \&i and ^ which should be interpreted as the boundary conditions. 
To get a nice expression for f'(<fi) from the Friedmann equation (|64p . we approximate <j){a) by the constant and 
the logarithmic term in (|66p . This is justified when w < 1, so that the second term decays fastest, and when a is 
large enough. Similarly, if we assume again large enough a, and the effective expansion is nonphantom when the 
physical matter content doesn't decay faster than dust, or more generally w > 2w rn — 1, then the solution (|67p may 
be approximated by setting \&2 = 0. Then we obtain 

3(w - w m )(3 -w- 2w m )<D? m a 2 +^ w ~ w ™) + ( w - 1)(1 + 3w)*i 
2(1 - 3w) 2 d 
and by integrating furthermore 

-2(3 -w- 2w m )n° m a 2 + 3 ^~^ +> - 1)(1 + 3^ 
4 (1 + w(3w — 4)) a z 

where the integration constant should be dropped as it corresponds to a constant shift of k 2 . By inverting the relation 
Eq. (|66|) we finally get the function as 

= I7TT k Tvi f( w - !)( 3w + l)*ie^ (0 - 0o) + 2(w + 2w m - 3)tl° m e* ■ (70) 

4(1 + w(6w — 4)) L J 

The main result about this example is that expansion with a constant equation of state (with or without matter) 
corresponds to an exponential coupling. Details are discussed below. 
We find two possibilities: a vacuum solution and a scaling solution. 

• In the case that w < w rn — 2/3 the asymptotic form of the matter coupling becomes tdominated by the first 
term inside the square brackets at large times. Thus, if we assume the effective equation of state to be more 
negative the ^-term dominates at large a instead of the matter term. Since the matter has diluted away, this 
corresponds to the vacuum solutions found by Jhingan et al (30j . As a function of a, the expansion rate of the 
vacuum corresponds to 

a — 1 

W e ff — . (71) 

With the addition of dust, these solutions become unstable. Thus one cannot reach them from the matter- 
dominated era. 

• If either the boundary condition is exactly = 0, or the expansion is less accelerating by w > w m — 2/3, 
then we find that the second term inside the square brackets dominates asymptotically and the solution is 
determined solely by the matter content. As one can deduce from Eq. (|57|) . in this solution the ratio Sl m /^ is a 
constant. In fact this is a scaling solution, since effective fractional matter density O m /(l + ^) is asymptotically 
a constant, since \1/ grows quickly to large values. This corresponds to a huge gravitational constant, which can 
be immediately excluded on observational grounds. The effective equation of state of this scaling solution is 
given as a function of a by 

1 , . 1 



w eff = 2^ + _ a _ g V 24 " - S6w m + (3 - 6a + 3w m ) 2 . (72) 

Since this is positive definite for w rn = and a < 0, one cannot use it to obtain accelerating from the dust 
dominated era. 
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10" 5 10" 3 10" 1 10" 5 10" 3 10" 1 

log(a) log(a) 

FIG. 1: Evolution of the scalar fields and comparison with the analytic approximations. The left panel is the power-law model 
with n — 3 (row 3 in Table |TJ) and the right panel is the exponential model with a — —0.5. The lower solid (blue) line is the 
numerical solution for (f> (divided by ten to fit into figure), and the dashed (black) line is the approximation (|4ip. The upper 
solid (cyan) line is the numerical solution for and the dashed (red) line the approximation (|5ip . We note that even after the 
universe has begun to accelerate, the matter-dominated solutions for the fields remain reasonable approximations. 




log(a) log(a) 

FIG. 2: The cosmological evolution leading to (nonaccelerating) effects from nonlocal corrections. 

Left panel: The power-law model with a = 6.0 and /o = 10 -6 . The solid (magenta) line is fi m /(l + *!/), and the dashed 
(red) line is the w e ff. One notes that the effect of / slowly decays and the universe asymptotically returns to the standard 
Einstein-de Sitter stage. For curiosity, we plot also the evolution in the Einstein frame. The dotted (black) line is Cl m , and the 
dash-dotted (blue) line is the w e f / , and These are understood as functions of the logarithm of the Einstein frame scale factor 
a. 

Right panel: The exponential model with a — —1.0 and fo — 10 -6 . The solid (magenta) line is fi m /(l + and the dashed 
(red) line is the w e ff. This is a scaling solution where we may consistently have fi m /(l + 9) > 1 because there can be negative 
effective energy sources. Here w e ff settles to the value given by (|72|) . The evolution is plotted also in the Einstein frame. The 
dotted (black) line is Cl m , and the dash-dotted (blue) line is the w e ff, and these are understood as functions of the logarithm 
of the Einstein frame scale factor a. In both frames we find asymptotic expansion with a nonzero constant effective equation 
state for the exponential model. 
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n 


fo 




4> 




1 


1.852216e - 02 


-9.701123e - 02 


-1.705136e + 01 


-6.316559e - 01 


2 


-1.068728e - 03 


-2.366542e - 01 


-1.673209e + 01 


-5.597827e - 01 


3 


6.090079e - 05 


-4.194428e - 01 


-1.656208e + 01 


-4.989559e - 01 


4 


-3.487378e - 06 


-6.412859e - 01 


-1.643555e + 01 


-4.461259e - 01 


5 


2.013315e - 07 


-8.990863e - 01 


-1.633576e + 01 


-4.010074e - 01 


6 


-1.172498e - 08 


-1.196474e + 00 


-1.625593e + 01 


-3.629566e - 01 


7 


6.875924e - 10 


-1.517736e + 00 


-1.618961e + 01 


-3.302487e - 01 


8 


-4.063656e - 11 


-1.891849e + 00 


-1.613562e + 01 


-3.027929e - 01 


9 


2.412112e - 12 


-2.252064e + 00 


-1.608829e + 01 


-2.783980e - 01 


10 


-1.440718e - 13 


-2.665776e + 00 


-1.604892e + 01 


-2.576894e - 01 



TABLE I: Numerical values for the power-law model with different exponents n. We report the scale fo which yields fi^ = 0.3. 
For even n this number must be negative. The following three columns then give the effective equation of the state and the 
values of the two fields evaluated today. 




log(a) log(a) 



FIG. 3: The cosmological evolution leading to the acceleration and singularity in various cases. 

Left panel: The power-law model with n = 3 and n = 6 (rows 3 and 6 in Table H} . The solid (magenta) line is fl m , and the 
dashed (red) line is the w e ff for n — 3. The dotted (black) line is Q m , and the dash-dotted (blue) line is the w e // for n = 6. 
The onset of acceleration leads to a singularity in the effective equation of state. The evolution is steeper for larger n. With 
smaller n, the relative amount of matter Q m , begins to decrease earlier and w e ff drops later to minus infinity. 
Right panel: The exponential model with a = —0.1 and a = —0.5. The solid (magenta) line is Q m , and the dashed (red) line 
is the Weff for a — 0.1. The dotted (black) line is Q m , and the dash-dotted (blue) line is the w e ff for a — 0.5. The onset of 
acceleration leads to a singularity in the effective equation of state. The evolution gets steeper for more negative a. To set 
tt m = 0.3 today fixes f = -0.296 for a = -0.1 and f Q = -0.140 for a = -0.5. 



IV. CONCLUSIONS 

We studied the dynamical evolution of cosmological expansion in a class of nonlocal gravity models extending 
general relativity by the addition of terms of the type /(i?/D), which are motivated by quantum corrections to the 
gravity action. The corrected action was rewritten as a multi-scalar-tensor theory, where several aspects of the theory 
such as its physical degrees of freedom, become more transparent. We also considered Weyl rescaling, but since the 
Einstein frame of these theories was not practical we performed our main analysis in the physical Jordan frame. The 
main result of this paper is the first explicit demonstration that these models can explain the present acceleration of 
the universe with simple forms of f(R/D) involving only Planck scale constants. 

In particular, power-law models and exponential models were used as examples. Schematically, we may summarize 
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the different possibilities of cosmological evolution in the exponential model as follows: 

/ = foe a * 

Matter domination Scaling, Eq. ([72]) Acceleration >■ Singularity 

Note that in the case a > one may have effects in the early universe. The power-law model features analogous 
possibilities. Now however the scaling is not exact, but the significance of the nonlocal correction decays and leads 
back to the usual Einstein - de Sitter stage. Schematically, we summarize this as: 

/ = M 

n<Q 

Matter domination Slows down expansion Acceleration >- Singularity 

Again, though this was not in the focus of the present study, n < may lead to a decaying effect which, with large 
enough scale /o, can influence the evolution of the early universe. Let us note that in both of the above cases, the 
conditions for acceleration to occur at late times imply f'(4>) > 0. According to Eq.©, this is precisely the condition 
for the field </> not to be a ghost. The consistency condition for the other field, W > — 1, ensuring the positivity of the 
graviton energy, is satisfied by all of the models we have considered. 

We sketched also a reconstruction method and applied it to find scaling solutions including matter and vacuum 
solutions with constant w e ff. These vacuum solutions have been found previously and shown to feature a possible 
de Sitter expansion [28| or effective dark energy with a range of equations of state extending to phantom regime [3C| . 
We confirmed and generalized these solutions, and noted that by addition of matter they become unstable, and thus 
cannot be reached in realistic cosmology. This shows that a detailed analysis of dynamics, without neglecting the effect 
of matter, is necessary in assessment of the cosmological relevance of these models. One could draw some analogy with 
previous findings about the cosmological evolution in the Gauss-Bonnet dark energy [43j], or in the gravity models 
of the f(R) type Q: though standard cosmological solutions could be (approximately) recovered at higher curvature 
and accelerating solutions would exist in vacuum at low curvature, the conditions for the cosmological dynamics to 
connect the first with the other may be very nontrivial [3, . To address the cosmological viability, it is also crucial 
to investigate the perturbations, as in the two classes of models mentioned above HESS]; for the present models this 
is left for future study. The form of the modifications is such that the impact to local gravity experiments could be 
expected to be negligible, but this should be checked in detail as well. 

A rather dramatic prediction of these models is that the present acceleration is a manifestation of our universe 
plunging into a sudden singularity. Though it awaits in the future and thus does not affect our present observations, 
its presence may be considered as an unpleasant property of the models signalling the need to choose more carefully the 
phcnomenological gravity action. In our tentative example, we showed explicitly the disappearance of the singularity 
by shifting □ — > D + eR. Finally, it is worth mentioning that the particular form of nonlocalities studied here is a direct 
generalization of an action previously considered in the quite different context of stabilizing the Euclidean gravity 0] . 
It is interesting that one may constrain nonperturbative quantum corrections to gravity by their cosmological effects 
at the present day, and that nonlocality may provide the required new physics to such seemingly detached major 
problems as the black hole information paradox and the fine-tuning of the cosmological constant. 



Nonlocal effect 



/o>0 




Nonlocal effect 



(-l)"/o>0 



(-i)"/o<o 
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